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ABSTRACT 

We investigate the relative occurrence of radio-loud and radio-quiet quasars in the 
first billion years of the Universe, powered by black holes heavier than one billion 
solar masses. We consider the sample of high-redshfit blazars detected in the hard 
X-ray band in the 3-years all sky survey performed by the Burst Alert Telescope 
(BAT) onboard the Swift satellite. All the black holes powering these blazars exceed 
a billion solar mass, with accretion luminosities close to the Eddington limit. For each 
blazar pointing at us, there must be hundreds of similar sources (having black holes 
of similar masses) pointing elsewhere. This puts constraints on the density of billion 
solar masses black holes at high redshift (z > 4), and on the relative importance of 
(jetted) radio-loud vs radio-quiet sources. We compare the expected number of high 
redshift radio-loud sources with the high luminosity radio-loud quasars detected in the 
Sloan Digital Sky Survey (SDSS), finding agreement up to z 3, but a serious deficit 
at ^ > 3 of SDSS radio-loud quasars with respect to the expectations. We suggest 
that the most likely explanations for this disagreement are: i) the ratio of blazar to 
misaligned radio-sources decreases by an order of magnitude above z = 3, possibly as 
a result of a decrease of the average bulk Lorentz factor; ii) the SDSS misses a large 
fraction of radio-loud sources at high redshifts, iii) the SDSS misses both radio-loud 
and radio-quiet quasars at high redshift, possibly because of obscuration or because 
of coUimation of the optical-UV continuum in systems accreting near Eddington. 
These explanations imply very different number density of heavy black holes at high 
redshifts, that we discuss in the framework of the current ideas about the relations of 
dark matter haloes at high redshifts and the black hole they host. 

Key words: BL Lacertae objects: general — quasars: general — radiation mecha- 
nisms: non-thermal — X-rays: general 



1 INTRODUCTION 

Ajello et al. (2009, hereafter A09) recently published the list 
of blazars detected in the all sky survey by Swift/B AT , be- 
tween March 2005 and March 2008. BAT is a coded mask 
telescope designed to detect Gamma Ray Bursts (GRBs), 
has a large field of view (120° x 90°, partially coded) and is 
sensitive in the [15-150 keV] energy range. This instrument 
was specifically designed to detect GRBs, but since GRBs 
are distributed isotropically in the sky, BAT performed, as a 
by product, an all sky survey with a reasonably uniform sky 
coverage, at a limiting sensitivity of the order of 1 mCrab in 
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the 15-55 keV range (equivalent to 1.27 x 10~^^ erg cm~^ 
s~^) in 1 Ms exposure, see A09). Taking the period March 
2005 - March 2008, and evaluating the image resulting from 
the superposition of all observations in this period, BAT 
detected 38 blazars (A09), of which 26 are Flat Spectrum 
Radio Quasars (FSRQs) and 12 are BL Lac objects, once 
the Galactic plane (\b < 15°|) is excluded from the analysis. 
A09 reported an average exposure of 4.3 Ms, and considered 
the [15-55 keV] energy range, to avoid background prob- 
lems at higher energies. The well defined sky coverage and 
sources selection criteria makes the list of the found blazars 
a complete, flux limited, sample, that enabled A09 to cal- 
culate the luminosity function (LF) and the possible cosmic 
evolutions of FSRQs and BL Lacs, together with their con- 
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tribution to the hard X-ray background. A09 also stressed 
the fact that the detected BAT blazars at high redshift are 
among the most powerful blazars and could be associated 
with powerful accreting systems. Within the BAT sample, 
there are 10 blazars (all FSRQs) at redshift greater than 
2, and 5 at redshift between 3 and 4. All (and only) these 
blazars have a X-ray luminosity exceeding Lx — 2 x 10*^ 
erg s~^. All these sources have been studied by Ghisellini et 
al. (2010, hereafter GIO), that showed that their optical-UV 
emission is dominated by the emission of their accretion disk, 
with no contamination from the beamed non-thermal con- 
tinuum, even if the latter is dominating the total bolometric 
luminosity. Fitting the optical-UV emission with a standard 
Shakura & Sunyaev (1973) accretion disk, it was possible to 
estimate both the mass and the disk luminosity. These high 
redshift blazars are shown in Fig. [1] (diamonds) together 
with all the blazars with z > 2 detected above 100 MeV 
in the 11 months all sky survey performed by Fermi/LAT 
(Abdo et al. 2010) and studied in Ghisellini et al. (2011). 
Fig. [1] shows that all high redshift BAT blazars are charac- 
terized by large Eddington ratio of order 0.2-1 and by black 
holes heavier than l(f Mq 

Since these objects are at high redshifts, our finding has 
important implications on the number density of heavy black 
holes, especially if we consider that for each blazar pointing 
at us, there must be hundreds of similar sources (having 
black holes of similar masses) pointing elsewhere. In fact, if 
the emitting plasma is moving with a bulk Lorentz factor 
r in one direction, the number of sources observed within 
the beaming angle l/F is only a fraction l/(2r'^) of the the 
sources pointing in other directions. 

Taking the luminosity function of A09 at face value and 
calculating the expected number of luminous sources (i.e. 
Lx > 2 X 10*^ erg s~^, likely hosting a black hole with 
M > Mq, where Mg = W^Mq) at z > 4, one finds that 
the number density of heavy black holes of jetted sources 
is close to or even greater than the upper limit defined by 
standard "dark matter halo-black holes" relationships at the 
largest redshifts. GIO then corrected the original A09 lumi- 
nosity function by assuming an evolutionary model that is 
equal to the A09 one up to 2 ~ 4.3 (where they measure the 
peak of the density of high X-ray luminosity blazars), and 
then cuts off exponentially ("minimal" LF). This "minimal 
LF" was consistent with the constraints posed by standard 
"dark matter halo-black holes" relationships and consistent 
with the constraints given by the existence of a few blazars 
(discovered serendipitously) at 2 > 4 (see below). 

In this paper we explore the implications of GIO results 
in view of the properties of the radio quiet and radio-loud 
populations, their redshift evolution and their connection to 
host dark matter halos. First we check if the expected num- 
ber of radio-loud sources calculated from the A09 and GIO 
luminosity functions agrees with those detected in the SDSS 
survey of quasars. As discussed in Section 2, while a rough 
agreement is found up to z ~ 3, there is a serious deficit 
of SDSS radio-loud sources above this redshift. We then in- 
vestigate the possible reasons of this discrepancy, discussing 
three possible solutions. While we cannot confidently select 
one of these, we point out the consequences they have on our 
understanding of the physics of jets and on the relationship 
between dark matter halos and the mass of the black holes 
they host. 
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Figure 1. Accretion disk luminosity Ld as a function of black 
hole mass for blazars with z > 2 in the BAT sample (diamonds; 
A09) and in the ILAC Fermi/LAT sample (circles, Abdo et al. 
2010), as studied in GIO and in Ghisellini et al. (2011). The high 
redshift BAT blazars are all characterized by black holes with 
M > IO^Mq and by Ld/^Edd > 0.1. 

2 RADIO-LOUD HIGH REDSHIFT SOURCES 

One can estimate the volume density of high redshift blazars 
hosting a black hole of mass larger than Mg using the cos- 
mological evolution model of A09 along with its high-z cut- 
off (i.e. "minimal") version, assuming that all blazars with 
Lx > 2x 10**'' erg s"* have & M > M9 black hole (GIO). 
We cannot exclude that blazars with lower X-ray luminosity 
also host massive black holes, so the "observational" points, 
strictly speaking, are lower limits. 

Lower limits to the density of high redshift blazars pow- 
ered by black holes with M > Mg are placed by the exis- 
tence of at least 4 blazars at 4 < 2: < 5 for which GIO 
have estimated a black hole mass larger than 10^ Mq. These 
blazars are RXJ 1028.6-0844 {z = 4.276; Yuan et al. 2005); 
GB 1508+5714 (z = 4.3; Hook et al. 1995); PMN J0525- 
3343 {z = 4.41; Worsley et al. 2004a) and GB 1428-H4217 
{z = 4.72; Worsley et al. 2004b). The lower limit in the 5-6 
redshift range corresponds to the existence of at least one 
blazar, Q0906-1-6930 at z = 5.47, with an estimated black 
hole mass of 2 X 10^ M0 (Romani 2006). 

These are all sources pointing at us. The real number 
density of heavy black holes in jetted sources, $rl(z, Af > 
Mg), must account for the much larger population of mis- 
aligned sources. We have then multiplied the mass function 
of blazars by 2r^ = 450, i.e. we have assumed an average 
F-factor of 15, appropriate for the BAT blazars analysed in 
GIO. 

Summarizing, the BAT blazar survey allowed to mean- 
ingfully construct the hard X-ray LF of blazars. GIO have 
also constructed the minimal evolution consistent with the 
existing data and the (few) existing lower limits. At the high 
luminosity end the LF can be translated into the mass func- 
tion of black holes with more than one billion solar masses. 
In Fig.[2]we show "I>r,l(2, M > Alg) as derived from the cos- 
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Table 1. " Redshift bin. ''Number of uniformly selected quasars with logLtoi > 47 (erg s~^) in the SDSS+FIRST survey (8770 square 
degrees in common). " Number of objects with the radio (5 GHz) to optical (2500 A) monochromatic flux ratio larger than 10. Column 
c divided by column b. " Number of radio loud objects expected in the SDSS+FIRST footprint, ~8770 square degrees. The number is 
obtained from the SDSS QSO luminosity function of Hopkins et al. (2007) rescaled by the radio loud fraction given in column d. ^ The 
number of blazars detected by Swift/3AT with an estimated black hole mass exceeding 10^ Mq, from the minimal and the A09 evolutions, 
when relevant. " The expected number of radio-loud sources calculated by multiplying the number of detected blazars (column /) in the 
given redshift bin by 450 X 8, 770/40, 000 ~ 10^ (i.e. assuming F = 15). Same as column g, but for F = 5 (the number scales as F^). 
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Figure 2. Number density of black holes with M > 10^ AIq as 
a function of redshift. The filled black square in the 3 < ^ < 4 
bin is taken directly from Fig. 10 of A09. The two empty squares 
account for the population of misaligned sources, multiplying by 
2F2, with F = 5 and F = 15. Red hatched area: number density of 
heavy black holes in radio-loud sources derived from the blazar 
LF studied in A09, assuming F = 5 (lower bound), or F = 15 
(upper bound). Green hatched area: number density of heavy 
black holes assuming the "minimal" cosmic evolution of the LF 
of blazars (GIO), with F = 5 (lower bound), or F = 15 (upper 
bound) . Blue hatched area: number density of heavy black holes in 
radio— quiet quasars (assuming the LF and its evolution of H07), 
assuming an average Eddington fraction of f^dd = 0-3 (upper 
boundary) or f^dd = 1 (lower boundary). Blue stars: number 
density of M > 10^ Mq derived from the existence of the black 
holes analyzed in Kurk et al. (2007). Purple dot: number density 
of M > 10^ Mq derived from the mass function of black holes at 
z = 6 proposed by Willott et al. (2010). 



mological evolution model of A09 as a red stripe, and that 
derived from the "minimal" LF as a green stripe. 

The BAT blazars described above can be compared to 
the radio-loud sources in the quasar catalog of the Sloan 
Digital Sky Survey (SDSS, Schneider et al. 2010) Data Re- 



lease Seven (DR7) that includes information on radio detec- 
tion in and the Faint Images of the Radio Sky at Twenty-cm 
survey (FIRST, Becker et al. 1995). The region of the sky 
covered by both surveys is ~8770 square degrees. We adopt 
the public catalog with quasar properties described in Shen 
et al. (2010), which includes quasars bolometric luminos- 
ity (using bolometric corrections derived from the compos- 
ite spectral energy distributions from Richards et al. 2006). 
The catalog also provides the radio flux density at rest- 
frame 6 cm and the optical flux density at rest-frame 2500 
A that can be used to calculate the radio-loudness. Follow- 
ing Jiang et al. (2007) we define a source radio-loud if it has 
radio to optical flux ratio, R, larger than 10. For a hand- 
ful of sources where optical quantities are not provided, we 
supplement the "raw" catalog by calculating the bolometric 
luminosity from the absolute i-band magnitude, assuming 
a bolometric correction of 2.5. This bolometric correction is 
derived by matching the average bolometric luminosity pro- 
vided by Shen et al. (2010) with the bolometric luminosity 
calculated from the absolute i-band magnitude for sources 
where the catalog lists both quantities. We then calculate 
the rest-frame optical flux from the luminosity, in order to 
derive an estimate of the radio-loudness. This "extended" 
catalog will be our reference. 

We select all sources that are in the FIRST+SDSS foot- 
print, have an optical bolometric luminosity > lO*'^ erg s~^ 
and R > 10. We also require the quasars to be selected 
uniformly using the flnal quasar target selection algorithm 
described in Richards et al. (2002). These amount to, e.g., 
21 radio-loud sources in the 4 < z < 5 redshift bin. We note 
that this number is not representative of a complete, volume 
limited sample (i.e., of the true luminosity function). To de- 
rive a simple estimate of the statistical incompleteness, we 
compare the number of quasars in Tab. [1] (column b) to the 
number predicted by the bolometric LF, which is derived 
from surveys that include the SDSS (Hopkins et al. 2007, 
hereafter H07), in the SDSS+FIRST area. For example, in 
the same redshift bins, approximately a factor ~ 10 more 
objects are expected from the LF proposed by II07. 

On the other hand, such "incompleteness" biases should 
not affect number ratios, such as the radio loud fraction 
(RLF) derived firom our FIRST+SDSS sample. As a matter 
of fact, our derived RLF (column d of Tab. [1} is completely 
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consistent with the RLF derived by Jiang et al. (2007) . This 
allows us to estimate the "expected number" of RL objects 
simply multiplying the values obtained from the H07 LF by 
the RLF we have derived from our sample. Such numbers are 
reported in column e of Tab.[TJ and should be compared to 
the expectations from detected BAT blazars (columns g and 
h). It is immediately clear that, at least at 2;>3, the expec- 
tations largely exceed what derived from the observed LF if 
r = 15, while r = 5 seems to be quantitatively consistent 
with data. For increasing redshift, the corrected fraction of 
RL objects (column e) is progressively lower, and its ratio 
with respect to columns g and h is decreasing (the figure in 
column e exceeds both columns g and h for z = 2 — 3, it is 
in between the values of columns g and h for 2: = 3 — 4, it is 
lower for 2 = 4 — 5 and 2 = 5 — 6. In the latter case, the cor- 
rected fraction of RL objects just coincides with the lower 
boundary of the range from column h, that adopts F = 5). 
We note that a discrepancy between the observed number 
of radio-loud quasars and theoretical predictions was first 
noted by Haiman et al. (2004) and confirmed by McGreer 
et al. (2009). 

Before discussing the possible nature of such discrep- 
ancy, it is useful to analyze the predictions concerning the 
fraction of radio-loud sources at high redshifts. 



2.1 Radio— loud vs radio— quiet quasars at high 
redshift 

We now estimate the ratio of radio-loud to radio-quiet 
quasars having black holes exceeding M9 as a function 
of redshift. We must therefore derive the number density 
$rq(z, M > Mg) of radio-quiet quasars hosting black holes 
of M > Mg in different redshift bins. Consider then the 
bolometric LF of radio-quiet sources (H07) and simply as- 
sume that quasars radiate at an average fraction /Edd of the 
Eddington limit so that: 
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To derive $rq(2:, M > M9) one then integrates the LF above 
the luminosity threshold, Lmin, corresponding to M = Mg. 

Estimates of $rq(2:,M > Mg) for radio-quiet quasars 
are shown in Fig. [2] for /Edd = 0.3 (Lmin ~ IQ^^Lq) and 
/gdd = 1 (Lmin = 3.4 X lO^^L©). Note that the lower the 
average /Edd, the lower the luminosity of a quasar that hosts 
a M > black hole is. Therefore, decreasing /Edd allows 
one to integrate the LF down to lower luminosities, thus 
increasing the overall $rq(z, M > Mg). However, if we were 
to assume, say, /Edd = 0.1, then the range of luminosities 
where we compare radio-loud and radio-quiet quasars would 
differ, as in blazars the (observed) accretion disk component 
is always more luminous than several x lO**^ erg s~^ (cf. Fig. 

m- 

Finally, if the SDSS misses quasars because of obscura- 
tion biasing optical selection (see below and Treister et al. 
2011), then this mass function is in reality a lower limit, as 
more active black holes might exist. 

We quantify the redshift evolution of the ratio of radio- 
loud vs radio-quiet sources in Fig. [2] and Fig. [S] We stress 
that up to z = 4, where we do see blazars, the cosmological 
evolution model, as derived by A09, is secure. Beyond z = 4 
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Figure 3. Ratio of the number density of radio-loud to radio- 
quiet quasars with black hole masses exceeding 10^ Mq. Blue and 
grey symbols refer to the ratio calculated for the SDSS sources 
(this work and Jiang et al. 2007 respectively), while green (lower 
set) and red (upper set) dots refer to the ratio of the radio-loud 
population inferred from the Fermi (at z = 1.5) and BAT blazars 
and the radio-quiet quasars as derived from the SDSS (H07). Red 
points have been calculated assuming the A09 cosmic evolution 
of blazars; green points correspond to the "minimal" cosmic evo- 
lution of blazars suggested in GIO. The radio-quiet population 
is estimated from the LF of radio-quiet quasars and its redshift 
evolution, from H07. Top panel: average Eddington fraction of 
radio-quiet quasars is /Edd = 0.3. Bottom panel: average Ed- 
dington fraction of radio-quiot quasars is /Edd = 1- Left panels: 
r = 15. Right panels: F = 5. 



it depends strongly on the assumed evolution. Since, how- 
ever, the "minimal" evolution provides a lower limit to the 
number of radio-loud systems, we can be assured that the 
radio-loud vs radio-quiet fraction remains at least close to 
constant, and near unity, up until 2 ~ 6. We find that the 
fraction of jetted sources increases from 2 = 3.5 to 2 = 4.5 
by roughly an order of magnitude. Fig. [2] also shows that 
for M > h'h and L > lO^L© (L>10'''' erg s"^) the number 
density of radio-loud quasars approaches and possibly pre- 
vails over that of radio quiet-quasars, if we take face value 
the extrapolation of the cosmic evolution suggested by A09. 

We further check our results via a comparison of the 
radio-loud fraction that we derive from the FIRST-f SDSS 
sample we uniformly selected. This radio-loud fraction is 
shown in Fig. |3] (blue squares). We compare our estimate 
with the results for quasars of similar luminosities of Jiang 
et al. 2007, shown as gray triangles (their Fig. 3, lower left 
panel). The point at 2 = 1.5 is based on the Fermi blazars 
analyzed in Ghisellini et al. 2011, in order to obtain a radio- 
loud fraction at low redshift and maximize the range where 
we can compare our results to Jiang et al. (2007). The agree- 
ment between our selection in the SDSS-f FIRST and Jiang's 
is excellent where the analyses overlap (2 ^ 4). A striking 
result we find is that, while at 2 ^ 2.5 the "parent popu- 
lation" of SDSS-t-FIRST radio loud quasars traces almost 
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perfectly the BAT blazars, assuming F = 15 and /Edd = 1, 
the two selections deviate at higher redshift. 

This analysis, which relies on a uniformly selected sam- 
ple, does not require a volume complete sample, as we are 
now dealing with fractions. Wo still find a dearth of radio- 
loud sources at high redshift. 

We stress once again that at z < 2.5 the blazar popu- 
lation, with r — 15, joined with the radio-quiet population 
with /Edd = 1, is in excellent agreement with SDSS/FIRST 
data (both our analysis and .Jiang ct al. 2007 analysis). At 
z = 3.5 the number density of blazars is derived from ob- 
served sources (no redshift extrapolation) , and our only as- 
sumption is the value of F. We are therefore confident that 
there must be either a transition in the astrophysical prop- 
erties of the population or a selection bias. Below we discuss 
some possibilities to explain the found discrepancy. 



2.2 Where Eire radio— loud high— redshift quasars? 

The problem wc face is simply that BAT blazars indicate 
that there must be many radio-loud sources that instead 
the SDSS-I-FIRST survey does not detect. In the following 
we list possible solutions. 

(i) The value F --^ 15 for the average bulk Lorentz factor 
of blazars is too large. A value of F ~ 5 would make the 
predicted numbers of misaligned radio sources to decrease by 
an order of magnitude (it is proportional to F^), becoming 
then consistent with the SDSS-I-FIRST detected radio-loud 
sources. We have checked that fitting the SED of our blazars 
with F = 5 gives reasonable results, but this value of F 
cannot account for all the measured apparent superluminal 
velocity and, furthermore, it implies that the jet is away 
from the conditions of minimum jet power requirements (see 
GhiscUini & Tavccchio 2010). Since the discrepancy appears 
as redshift increases, we need that F evolves with cosmic 
time (F must decrease with increasing redshift). 

A similar solution is that the jet has a radial velocity 
structure, similar to Gamma Ray Bursts, thus emitting (at 
approximately the same level) within an angle larger than 
1/F. We thus observe preferentially that part of the source 
exactly pointing at us (i.e. with a viewing angle close to zero) 
maximally beamed. The other slightly misaligned part of the 
source contribute less to the observed fiux. This implies that 
we underestimate the jet power, that refers to the part of 
the jet mostly contributing to the fiux, the one pointing at 
us. The other parts, viewed at angles larger than 1/F, are 
not accounted for when estimating the jet power. In other 
words, we calculate the power of only a part of the jet, of 
solid angle ^ 1/F'^, which is smaller than the jet solid angle 
Oj ~ df, where 9j is the jet half opening angle. Therefore, 
very approximately, we underestimate the jet power by the 
factor (Oj/r)'^. To account for the disagreement in number, 
(Sj/F)^ should be one order of magnitude. 

(ii) There is a bias, in the SDSS-l-FIRST survey, against 
the detection of powerful radio-loud sources at high-«, but 
not against radio-quiet quasars. This (yet unknown) bias 
could be due to the compactness of the radio halo in radio 
sources that are still too young to have developed an ex- 
tended structure. As a consequence, the fiux emitted in this 
compact and isotropic structure is self-absorbed up to the 



GHz frequency range. Further study is however needed to 
verify this possibility, that now is only a speculation. 

(iii) The SDSS/FIRST selection misses a large fraction 
of powerful quasars (both radio-loud and radio-quiet) at 
z > 3. This may be the results of optical absorption, or 
else of coUimation of the optical emission of the disk (mak- 
ing the apparent disk luminosity much dimmer if the disk 
is observed edge on). An absorbed (or optically anisotropic) 
source intrinsically emitting more than 10*'^ erg s^^ in the 
optical (our selection limit) would not pass our selection. 
This would occur also for radio-loud objects, because, even 
if the radio-emission is unaffected, the primary selection 
(SDSS) is on the optical luminosity. 

The radio-loud fraction could be right (if the bias applies 
equally to both kind of sources), but both classes are under- 
represented by an order of magnitude, at least at « > 3. 

These possibilities arc listed in order of increasing de- 
mands of heavy black holes at large redshifts. In fact, case i) 
would minimize the total number of high-z powerful radio- 
loud sources (and their associated heavy black hole), that 
would be well described by the SDSS-I-FIRST survey. The 
price to pay is a correspondingly larger requirement on the 
jet power of each quasar (i.e. we have a factor 10 less radio- 
loud powerful sources, and then less heavy black holes, but 
each jet is 10 times more powerful). In this case the total 
number of heavy black holes at high-z is practically given 
by the radio-quiet quasars (the radio-loud contributing by 
less than 10%). 

Case ii) is intermediate, since it implies that the num- 
ber of powerful radio-loud sources is correctly described by 
the BAT blazars (with (F) ~ 15), that should be as numer- 
ous as the radio-quiet sources (with M > A/g) detected by 
the SDSS. Therefore the number of heavy and early black 
holes is roughly twice as much as the one derived by the 
radio-quiet quasars LF. The most important consequence 
in this case is the strong evolution of the radio-loud frac- 
tion, becoming close to unity at 2 > 3, with interesting con- 
sequences on our understanding of the growth of the early 
supermassive black holes. 

Case iii) is the more demanding: it implies that the 
SDSS misses 90% of the most powerful quasars (no matter 
their radio loudness) and thus that the number of heavy 
black holes is a whole order of magnitude more than the 
mass function derived from the H07 LF would predict. 

In the following we will examine the current ideas about 
the relationship between dark matter halos and black hole 
mass, in order to estimate the expected number of heavy 
black holes as a function of redshift. 



3 BLACK HOLE-DARK HALO CONNECTION 

Empirical correlations have been found between the central 
stellar velocity dispersion and the asymptotic circular veloc- 
ity (K) of galaxies (Ferrarese 2002, Baes et al. 2003, Pizzella 
et al. 2005): 

/ V \ ^-^^ 

o- = 200kms~M — — (2) 

V320kms-i J ^ ' 
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and 



: 200 km s 



339 kms~ 



(3) 



as suggested by Pizzella et al. (2005) and Baes et al. (2003), 
respectively. Some of these relationships (Ferrarese 2002, 
Baes et al. 2003) mimic closely the simple a — Vc//3 defi- 
nition that one derives assuming complete orbital isotropy. 
We note that in an isothermal sphere a = Vc//2. 

Since the asymptotic circular velocity (K) of galaxies 
is a measure of the total mass of the dark matter halo of 
the host galaxies, one can relate in simple ways the mass of 
the central black hole to the mass of its host halo ("hole- 
halo" connection, e.g., Ferrarese 2002, Volonteri et al. 2003, 
Wyithe & Loeb 2003, Rhook & Wyithe 2005, Croton 2009). 
A halo of mass Mh collapsing at redshift z has a circular 
velocity: 
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where Ac is the over-density at virialization relative to the 
critical density. For a WMAP5 cosmology we adopt here the 
fitting formula Ac = IStt^ + 82d - 39d^ (Brian & Norman 
1998), where d = Sl^ — 1 is evaluated at the collapse redshift. 
In this case we obtain 
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We will further assume that the black hole-cr (M-ct) scaling 
is: 

M _ / cr 

Mg ^ V356 kms- 

(Tremaine et al. 2002). We also assume that these scaling 
relations observed in the local universe hold at all redshifts. 
Therefore we derive relationships between black hole and 
dark matter halo mass: 
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if we adopt the relationship in Pizzella et al. (2005); 
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if we adopt a = Vc//3 (almost equivalent to what one would 
calculate using the relationship in Baes et al. 2003). 
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Figure 4. Mass functions of black holes at 2 = (top) and at 
2 = 6 (bottom). Vertical (red) bar: constraints at 2 = (Giiltekin 
et al. 2009). Dash-dotted (blue) curve: a = Vc//3 (Eq. 8)+ M-cr; 
short dashed (orange) curve: Pizzella (Eq. 7) +M-a. As a refer- 
ence, we also show also the mass function obtained by assum- 
ing M = /bar lO^^Mh (black long— dashed curve) and the case 
M = 10~^Mii (green dotted curve). 
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3.1 Black hole mass functions 

We can therefore estimate the mass function of black holes 
by convolving Eq. 7 and Eq. 8 with the mass density of 
dark matter halos with mass Mh derived from the Press & 
Schechter formalism (Sheth et al. 1999). The number den- 
sity of black holes with M > Mg, (^{z,M > Afg), there- 
fore corresponds to the number density of halos with mass 
Mh > Mthr, if Mthr is the mass of a halo that hosts a billion 
solar masses black hole, and we assume that all halos host 
black holes. 

In Fig. |4] we show the mass functions derived via Eq. 
7 and Eq. 8 coupled to the Press & Schechter function at 



Figure 5. Number density of black holes with M > 10® Mq as 
a function of redshift. Red hatched area: number density of black 
holes in radio-loud sources, derived from the blazar LP of A09, 
assuming F = 5 (lower bound), or F = 15 (upper bound). Green 
hatched area: "minimal" number density (studied in GIO), assum- 
ing with F = 5 (lower bound), or P = 15 (upper bound). Blue 
hatched area: black hole number density of radio-quiet quasars 
(from the LF and its evolution studied in H07) Line styles and 
colors of the "hole-halo" connection predictions as in Fig. |4] 



2 = (bottom panel) and z — 6 (top panel). As an exer- 
cise, one can also assume that the black hole mass scales 
linearly with the halo mass. We can derive a plausible up- 
per limit to this scaling assuming M = /blO"^Mh, where 
/b ~ rib/^M — 0.14 is the universal baryon fraction, and 
M ~ 10~'^Mbuige is the empirical correlation between black 
hole and bulge mass in elliptical galaxies (Marconi & Hunt 
2003, Haring & Rix 2004). This assumption corresponds to 
assuming that galaxies do not loose any baryon because of 
feedback effects, and all baryons end up in a stellar bulge. 
At z = this relationship is obviously wrong (the baryon 
content in galaxies is much less than /b and not all baryons 
end up in stellar bulges), and we find that M — 10~^Mh 
provides a more acceptable solution, as shown in the upper 
panel of Figure ID At 2: = the mass function of black holes 
is estimated in the literature by coupling the empirical cor- 
relations found between black hole mass and host properties 
(bulge mass, luminosity and velocity dispersion, see Marconi 
et al. 2004, Gultekin et al. 2009 and references therein) with 
the distribution of galaxies as a function of these proper- 
ties. In the top panel of figure [4] we show with a vertical bar 
the current limits on the mass density of black holes with 
M > Mg at 2; = 0. The "hole-halo" connection coupled to 
the Press & Schechter function might possibly be slightly 
overestimating the number density of large black holes at 
z = 0. This is due to the LF of galaxies being steeper at the 
high mass end than the halo mass function at the high mass 
end. 

We stress here that at least two biases do exist that 
affect the comparison between the "hole-halo" connection 
and observational samples: first, Lauer et al. (2007) suggest 
that the M — a relation compared to M — L relation under- 
estimate the number of BH with M > Mg (but see Bernardi 
et al. 2007; Tundo et al. 2007). Second, going in the same 
direction, the intrinsic scatter in the M — a relation allows a 
larger number of possible haloes hosting massive black holes 
(see Lauer et al. 2007b; Gultekin et al. 2009). We discuss 
the importance of the intrinsic scatter in ^3.31 

Finally, Kormendy et al. (2011a) question any corre- 
lation between black holes and dark matter halos (but see 
Volonteri et al. 2011). We notice that Kormendy's argument 
is not a concern here, as at large masses Kormendy et al 
(2011b) suggest that a 'cosmic conspiracy' causes a and Vc 
to correlate, thus making the link between M and Vc ad- 
equate. Although estimates we derive from the halo-hole 
connection are therefore extremely uncertain, they can still 
provide some sense of the possible hosts of these massive 
black holes at high redshift. 

3.2 Number density of high-redshift M > W^Mq 
black holes powering blazars 

We now turn to compare the number density of M > Mg and 
L > 10"Lo (i>10'''^ erg s"^) blazars to the upper limits 
defined by "hole-halo" connection. The various assumptions 
for the "hole-halo" connection discussed in §3 are shown in 
Fig. [5l The number density of heavy black holes powering 
jetted sources is now close to or even greater (if we extend 
the cosmological evolution model of A09 beyond z~4) than 
the upper limit defined by "hole- halo" connections at the 
largest redshifts. The mass function derived by the "mini- 
mal" LF is instead consistent. 
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Figure 6. Importance of scatter for the prediction of the number 
density of black holes with M > 10^ Mq as a function of redshift. 
Top panel: the a = Vc/ZS (Eq. 8)-!-^ — a curve (blue dashed 
dotted) increases by assuming a scatter of 0.3 and 0.5 dex (dashed 
blue curves). Bottom panel: even the M = 10~^Mii curve (green 
dotted) can become consistent with the LF of radio-loud sources 
(hatched areas) assuming a scatter of 0.5 dex. The grey shaded 
area in both panels indicate the limits on today's number density 
of M > Mg black holes. 

Despite the uncertainties, we can make some simple 
inferences on the "hole-halo" connection at high-redshift: 
the number density of M > Mg black holes cannot be be- 
low the limit imposed by the number density of blazars de- 
rived from the "minimal" evolution. For instance, it can be 
ruled out that at high redshift the black hole mass scales as 
M = 10~^Afh, if scatter is negligible. This implies that high 
redshift black holes represent a higher fraction of the mass 
of a dark matter halo (at least, this is the case for the most 
massive active black holes, M > Mg and L > 10^^ Lq, this 
might not be the case for lower mass holes, see, e.g. the dis- 
cussion in Willott et al. 2010). Including the effect of scatter 
eases such constraints, as we show below. 

3.3 Importance of scatter 

We can assume that at fixed a the logarithmic scatter in 
black hole mass is A =0.3-0.5 dex (Mbh = Mbh,,t x 10"^^ 
where 5 is normally distributed, see, e.g., Gultekin et al. 
2009). We include scatter, at various levels of A, by perform- 
ing a Monte Carlo simulation, where for each black hole mass 
we create 500 realizations of the host mass. Fig. [6] shows 
two examples of Monte Carlo realizations of the mass func- 
tion and number density that include scatter at the level of 
A =0.3-0.5 dex (Merloni et al. 2010). By comparing Fig. [6] 
to Fig. [5] one clearly sees how scatter dramatically increases 
the number of black holes with M > Mg, and for most 
"hole-halo" connection the number density of high-redshift 
radio-loud sources can be accommodated. For instance, the 
M — a relations can accommodate the radio-loud population 
as long as scatter is around 0.3 dex. Notice that in the log- 
arithmic scale of Fig. |6] today's number density of M > Mg 
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number of massive black holes that have to exist at high 
redshift, as we have to take into account radio-quiet sources, 
for instance including optically selected quasars as described 
in Section 2.1. If the bolometric LF of radio-quiet sources 
(H07) is a good tracer of the quasar population, than there 
are roughly 0.1 radio- loud quasars per each radio-quiet one 
(assuming the minimal evolution of the LF of blazars), up 
to one or more, if we take the evolution of A09 face value. 

However, based on the arguments of §2.1 and §2.2, there 
might be a large population of radio-quiet quasars that 
are not accounted for in the bolometric LF of radio-quiet 
sources of H07 (see case iii in §2.2). Since we have no infor- 
mation on this putative hidden population we do not here 
include this speculation in our analysis, but if in reality most 
quasars are missed by current surveys, then all problems we 
discuss below are exacerbated. 

We can assess the requirements on the population of 
high-redshift BH hosts by investigating the fraction of M > 
Mg black holes that are actively accreting, that is, the active 
fraction of black holes. We define the active fraction as: 



Figure 7. Fraction of M > 10^ Mq that is active. Red: A09 cos- 
mic evolution inferred from the LF of blazars (this curve is basi- 
cally unaffected by different choices of /Edd '^^ radio-loud quasars 
dominate the active population). Green: "minimal" cosmic evo- 
lution inferred from the LF of blazars. Stars: radio-loud systems 
only. Filled circles: all active quasars, where we assumed /Edd = 1- 
Empty circles: all active quasars, where we assumed /Edd = 0.3. 
Bottom: a = VcA/3 + M -a. Top: Pizzella + M - ct. In all cases 
we ignore scatter in the M ~ a relation. The grey shaded area is 
not permitted as the active fraction becomes larger than unity. 



black holes (red vertical bar in the upper panel of Fig. 4) 
is between 3 and 5 (grey stripe in Fig. |6}, comparable to 
the number density at ^ ~ 5 in the cases with significant 
scatter .X For instance, if we assume the scaling of Eq. 8, the 
number density of M > Mg black holes reaches 10'' Gpc""' 
at 2 = 4 for A =0.3 and 2 = 7 for A =0.5. If A =0.5 the up- 
per limit to today's number density, 10^ Gpc""', is reached 
at 2 = 2, implying that after that cosmic time the number 
density of M > A/g black holes cannot increase any more. 

We have hitherto assumed that the duty cycle, a;dc = 1, 
corresponding to the fraction of black holes that are active 
(related to the ratio of the lifetime of quasars to the Hub- 
ble time), that is, all halos host an active black hole. We 
have further assumed that radio-loud quasars dominate at 
high redshift. If the active fraction a^dc is less than unity, im- 
plying that not all black holes with M > Mg are active and 
accreting close to the Eddington rate, than the requirements 
become stricter. Similarly, if radio-quiet quasars dominate 
the population of active black holes, their number has to 
be accounted for as well. In the next section we expand our 
analysis to include radio-quiet sources and allow for an ac- 
tive fraction, or duty cycle, below unity. 



3.4 Active fraction of high redshift M > 10^ M© 
black holes 

Until now we have focused only on the constraints that high- 
redshift blazars impose on the number density of Al > Mg 
black holes. However, this is clearly a lower limit to the 
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for radio-loud blazars only, and 



Xdc(RL-l-RQ) 



'1'rl(2, M > Mg) + $rq(2, M > Mg) 
^{z,M > Mg) 



(10) 

for all active sources, where $r,l(2,M > Alg) and 
'1>rq(2, M > Mg) are defined in section 2. <E>rq(2, M > Mg) 
depends on one parameter, the typical Eddington ratio of 
radio-quiet black holes. $(2,M > Mg) is defined in §3.1, 
and depends on the "hole-halo" connection and on the level 
of scatter that this relationship suffers. The numerator and 
denominator in both Eq. 9 and 10 are derived independently, 
hence a priori the active fraction can apparently assume val- 
ues above unity. When Xdc > 1 the result is however unphys- 
ical, and it allows us to rule out a given model. 

The active fraction is shown in Fig. [7] for /Edd = 
1, /Edd = 0.3 and negligible scatter in the "hole-halo" 
connection. We here show both the total active fraction, 
2;dc(RL-|-RQ), and the active fraction in radio-loud sources 
only, a;dc(RL). In both cases we have assumed F = 15. Face 
value, by 2 ~ 5 almost all "hole-halo" connections, except 
for the case of Eq. 7, require an active fraction (or a duty 
cycle) of unity. A significant amount of scatter might how- 
ever alleviate the issue, as shown in Fig. [8] As discussed 
above, scatter increases very significantly the number den- 
sity in M > Mg black holes, and pushes the active fraction 
to much lower values, at the cost however of having already 
built up almost all M > Mg black holes by 2 ~ 5 (cf. Fig. 6). 

Finally, if black holes accreted at rates significantly be- 
low the Eddington rate, the estimate of the number density 
of radio-quiet quasars that we derive from the luminosity 
function would increase (§3). This decreases the radio- loud 
fraction (see Fig. 5), but at the same time the total active 
fraction increases, and the active fraction becomes close to 
unity even for a significant scatter of 0.3 dex. 
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Figure 8. Same as Fig.Q for the case cr = Vc/ZS + M - 
including scatter. 
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4 DISCUSSION AND CONCLUSIONS 

We investigated the relative occurrence of radio-loud and 
radio-quiet quasars in the first billion years of the Universe, 
based on the sample of high-redshift blazars detected in the 
3-years all sky survey performed by the Burst Alert Tele- 
scope (BAT) onboard the Sm/t satellite (Ajello et al. 2009). 
The masses M of the black holes powering these quasars 
exceeds a few billions solar masses, with accretion luminosi- 
ties being a large fraction of the Eddington limit (GIG). For 
each blazar pointing at us, there must be hundreds of sim- 
ilar sources (having black holes of similar masses) pointing 
elsewhere. This can set constraints on the number density 
of dark matter halos that can host massive black holes at 
high redshift. 

We first compared the number of radio-sources hosting 
heavy black holes estimated from the BAT detected blazars 
to the SDSS4-FIRST survey to explore the relative impor- 
tance of (jetted) radio-loud vs radio-quiet sources. We find 
a rough agreement up to z ~ 3, but beyond this redshift 
there is a deficit of radio sources (detected by the FIRST 
and present in the SDSS surveys) with respect to the expec- 
tations (see also Haiman et al. 2004, McGreer et al. 2009). 
We found no obvious explanation for this deficit, and have 
suggested three possibilities: i) the bulk Lorentz factor of 
the jet (controlling the number of misaligned sources) de- 
creases beyond z = 2> (from F ~ 15 to F ~ 5); ii) there is a 
bias against detecting distant (and therefore possibly young) 
radio-sources with the FIRST survey (i.e. at 1.4 GHz) and 
iii) there is a bias against optical selection of distant and 
powerful quasars, both radio-quiet and radio-loud, due to 
absorption or coUimation of the disk emission. These possi- 
biUties aflect our estimates of the number density of heavy 
black holes (M > 10^ M0) in an increasing way [from possi- 
bility i) to iii)]. 

In the first case the majority of quasars are radio-quiet 
at all redshifts, and the number density of high-redshift 



M > 10^ Mq black holes can be safely derived from the ob- 
served LF of radio-quiet quasars. In the second case radio- 
loud and radio-quiet quasars powered by heavy black holes 
become comparable in number beyond z ~3-4, doubling the 
number of heavy and high-z black holes estimated from 
radio-quiet quasars only. This would also mean that the 
radio-loud fraction increases with z for sources with heavy 
black holes, suggesting that that a radio phase is perhaps a 
necessary ingredient for fast black hole growth at early cos- 
mic times. The last possibility is the most demanding, since 
it implies that we see only a minor fraction of the intrinsi- 
cally luminous high-z quasars (both radio-loud and quiet), 
implying that although the radio-loud fraction is always of 
the order of 10% (i.e. at all redshifts), the number of heavy 
black holes is now severely under-estimated (by one order 
of magnitude). 

We re-iterate that there is a good agreement between 
the number density of M > IO^A/q black holes found with 
blazars and the total number of radio-loud quasars up to 
2: ~ 3, but not beyond. We then conclude that even if we 
do not know the cause for this disagreement, there must be 
some change occurring at 2; ~ 3. 

We then studied plausible "hole-halo" connections in 
order to predict the number of supermassive black hole at 
high redshifts. We found that, unfortunately, the predicting 
power of these relations is weak, mainly because of the large 
effect that scatter can have: since we are dealing with fastly 
declining functions (corresponding to the high end of the 
distributions of luminosities and/or black hole masses), it 
is possible that even a few large black hole inhabiting halo 
slightly less massive than implied by the adopted relation 
can dominate the number density at a given redshift. 

On the other hand, despite the rather large uncertain- 
ties, the sources that have been already observed (cf. § 2) 
suggest that large and distant black holes are all active (or 
nearly so), and that they are all Eddington limited (or nearly 
so). If not, the number of heavy black holes would be larger, 
and the simple theoretical ideas we have adopted in this 
paper would start to have some difficulties to account for 
them. 

Finding more luminous radio-quiet quasars at z > 5 
is obviously mandatory to study the high mass end of the 
black hole mass density, but we would like to stress that 
finding high redshift blazars might be, in the end, even more 
important, since each one of those implies the existence of 
many more misaligned sources. A few blazars detected at 
z ~ 6 would be very challenging for structure formation, 
very constraining, and possibly illuminating for the under- 
standing the early growth of very massive black holes, and 
its feedback on the host. The existence of these blazars pos- 
sibly implies that normal "feedback" might not be at play 
at the highest redshifts. A possible explanation is that high- 
accretion rate events, distinctively possible during the vio- 
lent early cosmic times, trigger the formation of coUimated 
outfiows that do not cause feedback directly on the host, 
which is pierced through. These jets will instead deposit 
their kinetic energy at large distances, leaving the host un- 
scathed. This is likely if at large accretion rates photon trap- 
ping decreases the disk luminosity, while concurrently the 
presence of a jet helps dissipating angular momentum, thus 
promoting efficient accretion. This picture may explain why 
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high-redshift massive black holes can accrete at very high 
rates without triggering self-regulation mechanisms. 



ACKNOWLEDGMENTS 

We thank M. J. Rees for enlightening discussions. MV ac- 
knowledges support from SAO Awards TM9-0006X, TMl- 
12007X and NASA award ATP NNX10AC84G. RDC ac- 
knowledge financial support from ASI (grant n. 1/088/06/0). 



REFERENCES 

Abdo A.A., Ackermann M., Ajello M. et al., 2010, ApJ, 
715, 429 

Ajello M., Costamante L., Sambruna R.M., et al., 2009, 

ApJ, 699, 603 (A09) 
Baes M., Buyle R, Hau G. K. T., Dejonghe H., 2003, MN- 

RAS, 341, L44 
Becker, R. H., White, R. L., & Helfand, D. J. 1995, ApJ, 

450, 559 

Bernard!, M., Sheth, R. K., Tundo, E., & Hyde, J. B. 2007, 

ApJ, 660, 267 
Bryan G. L., Norman M. L., 1998, ApJ, 495, 80 
Croton D. J., 2009, MNRAS, 394, 1109 
Ferrarese L., 2002, ApJ, 578, 90 

Ghisellini G., Delia Ceca R., Volonteri M. et al., 2010, MN- 
RAS, 405, 387 (GIO) 
Ghisellini G., Tavecchio F., 2010, MNRAS, 409, L79 
Ghisellini G., Tagliaferri G., Foschini L. et al., 2011, MN- 
RAS, 411, 901 

Giiltekin K., Richstone D. O., Gebhardt K., Lauer T. R., 
Tremaine S., AUer M. C., Bender R., Dressier A., Faber 
S. M., Fihppenko A. V., Green R., Ho L. C., Kormendy 
J., Magorrian J., Pinkney J., Siopis C., 2009, ApJ, 698, 
198 

Jiang, L., Fan, X., Ivezic, Z., Richards, G. T., Schneider, 
D. P., Strauss, M. A., & Kelly, B. C. 2007, ApJ, 656, 680 
Kormendy, J. & Bender, R. 2011, Nat, 469, 377 
Kormendy, J., Bender, R., & Cornell, M. E. 2011, Nat, 469, 
374 

Kurk, J. D., et al. 2007, ApJ, 669, 32 

Haiman, Z., Quataert, E., & Bower, G. C. 2004, ApJ, 612, 
698 

Haring N., Rix H.-W., 2004, ApJ, 604, L89 

Hook I.M., McMahon R.G., Patnaik A.R., Browne I.W.A., 

Wilkinson P.N., Iwrin M.J. & Hazard C., 1995, MNRAS, 

273, L63 

Hopkins P.F., Hernquist L., Cox, T.J., Robertson B. & 

Krause E., 2007, ApJ, 669, 67 (A07) 
Lauer T. R., Faber S. M., Richstone D., Gebhardt K., 

Tremaine S., Postman M., Dressier A., AUer M. C, Filip- 

penko A. V., Green R., Ho L. C, Kormendy J., Magorrian 

J., Pinkney J., 2007a, ApJ, 662, 808 
Lauer T. R., Tremaine S., Richstone D., Faber S. M., 2007b 

ApJ, 670, 249 
Marconi A., Hunt L. K., 2003, ApJ, 589, L21 
Marconi A., Risaliti G., Gilh R., Hunt L. K., Maiohno R., 

Salvati M., 2004, MNRAS, 351, 169 
McGreer, I. D., Helfand, D. J., & White, R. L. 2009, AJ, 

138, 1925 



Merloni A. et. al., 2010, ApJ, 708, 137 

Pizzella A., Corsini E. M., Dalla Bonta E., Sarzi M., Coc- 

cato L., Bertola F., 2005, ApJ, 631, 785 
Richards, G. T., et al. 2002, AJ, 123, 2945 
Richards, G. T., et al. 2006, ApJS, 166, 470 
Rhook K. J., Wyithe J. S. B., 2005, MNRAS, 361, 1145 
Romani R.W., 2006, AJ, 132, 1959 
Shakura N.I. & Sunyaev R.A., 1973 A&A, 24, 337 
Schneider, D. P., et al. 2010, AJ, 139, 2360 
Shen, Y., et al. 2010, larXiv:1006.5178l 
Sheth R.K. & Tormen G., 1999, MNRAS, 308, 119 
Tremaine S. et al., 2002, ApJ, 574, 740 
Tundo, E., Bernardi, M., Hyde, J. B., Sheth, R. K., & 

Pizzella, A. 2007, ApJ, 663, 53 
Yuan W., Fabian A.C., Celotti A. & McMahon R.G., 2005, 

MNRAS, 358, 432 
Volonteri M., Haardt F., Madau P., 2003, ApJ, 582, 559 
Volonteri, M., Natarajan, P., & Gultekin, K. 2011, 

■arXiv:1103.1644 

Willott C. J., Albert L., Arzoumanian D., Bergeron J., 

Crampton D., Delorme P., Hutchings J. B., Omont A., 

Reyle C, Schade D., 2010, AJ, 140, 546 
Worsley M.A., Fabian A.C., Turner A.K., Celotti A. & Iwa- 

sawa K., 2004a, MNRAS, 350, 207 
Worsley M.A., Fabian A.C., Celotti A. & Iwasawa K., 

2004b, MNRAS, 350, L67 
Wyithe, J. S. B., & Loeb, A. 2003, ApJ, 595, 614 



